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Abstract: We present a physically inspired generalization of equilib- 
rium response formulas, the fluctuation-dissipation theorem, to Markov 
jump processes possibly describing interacting particle systems out- 
of-equilibrium, following the recent work of [H [2]. Here, the time- 
dependent perturbation adding a potential V with small amplitude ht 
changes the rates W{x,y) for the transition x ^ y into 

Wt{x, y) = W{x, y) e^' i^yiy)-^^-)) 

as first considered by Diezemann, [1]; a, 6 are constants. We observe 
that the linear response relation shows a reciprocity symmetry in the 
nonequilibrium stationary regime and we interpret the connection with 
dynamical fluctuation theory. 



Dedicated to the 15th birthday of 
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1. Hurrah 

The present paper is devoted to an important theme of statistical 
thermodynamics and system theory. It deals with the response of a 
system to an external stimulus. More specifically, we are interested in 
the linear response to an energy impulse applied to nonequilibrium sys- 
tems and generalizing the celebrated fluctuation-dissipation theorem. 
That question comes up in the general construction of nonequilibrium 
statistical mechanics, but has possible applications in a variety of con- 
texts. It is not clear yet whether the answer that we give here and that 
was presented in a more restricted sense in the physics literature [H [2] , 
is operationally useful. We do however attempt some interpretation di- 
rected towards dynamical fluctuation theory. The mathematical origin 
of all this is the theory of stochastic processes, here in its most sim- 
ple representation for Markov jump processes on finite alphabets. The 
very fact that physically interesting relations can possibly be suggested 
already from the elementary mathematical theory of Markov processes 
is rather encouraging for the Markov field. That is then our contri- 
bution to the celebration of a young journal devoted to that subject, 
hurrah and many years to come. 

2. Response in equilibrium 

Relations between fluctuations, response behavior and dissipation in 
equilibrium systems have been obtained and applied throughout the 
development of statistical mechanics in the 20th century, [7j. Quite 
often textbooks treat linear response in a quantum mechanical context 
when applied to discrete systems such as spins or particles hopping on 
a lattice. The method for that equilibrium fluctuation-dissipation theo- 
rem is first order perturbation theory on time-dependent Liouville-von 
Neumann equations. 

Here we consider stochastic evolutions, Markov jump processes; we give 
a more probabilistic treatment and corresponding statistical mechani- 
cal interpretation. The mathematics remains elementary. 

2.1. Equilibrium dynamics. Let us consider a simple situation, which 
falls in the context of the present discussion. Take an Ising spin system 
on a finite graph (A, ~); at each vertex z G A there is a spin a{i) = ±1. 
A spin flip Markov dynamics on the configurations a G {+1, —1}"^ = K 
has possible transitions a — where (J-'ii) = cr{i) for j i and 
= is the new configuration with the spin flipped at vertex 

j. Physically, we imagine that there is a thermal reservoir perhaps in 
the form of lattice vibrations or of electronic degrees of freedom at- 
tached to the system so that for each transition a — > there is an en- 
ergy exchange U{a^) — U (a) with and an entropy flux [U (cr) — U {a^)]/T 
in the reservoir at equilibrium temperature T. There is no need here 
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to specify that energy function U{a). As long as we assume the system 
is only in contact with the outside world through this one tempera- 
ture bath, there should be a reversible stationary distribution p giving 
probabilities 

Pi-) = ^e-^""""' (1) 

to the spin configurations. We call it the equilibrium distribution. The 
reversibility is expressed by the condition of detailed balance 

W{a, a^) p{a) = W[a\ a) p{a^), for all j e A,a e K 

for the rates W{a, a^) of the transitions a ^ . Remark that all static 
properties of the system depend only on that part of the rates that is 
determined by the detailed balance condition. In other words, as long 
as the rates satisfy 

they will simulate the correct physical properties of the system. Ex- 
plicitly, the transition rates of the stochastic Ising model are 

W{a,a^) = ij{a,j) exp-^[U{a^) - U{a)] (2) 

and all produce the same equilibrium ([T]) independent of the prefac- 
tor il){a,j) = Tpi^a^ ,i) as long as it indeed does not depend on 
Moreover there is then a time-reversal symmetry for the stationary 
process: denoting by Pp the stationary Markov (equilibrium) process 
with stationary law p, we have equilibrium correlations 

which are functions of |t — s|. The brackets (■)cq denote an average over 
the equilibrium ensemble over all possible realizations of the stochastic 
process determined by ([I])-®. 

2.2. Perturbation and response. Suppose now that we start in 
equilibrium p at time t = but thereafter we slightly modify the dy- 
namics in a time-dependent way. For times s G [0,t] we switch on 
a magnetic field of small amplitude hs- That is the external stim- 
ulus by which we change the energy function U into U — hsV for 
V{a) = XliGA^*- How will the equilibrium system respond at time 
t > s, and does the choice of ip in the rates ([2]) make a difference? 
We look at the linear response 

= (g(t))eq + r d. KR'^^yit, S) + 0{h) 

Jo 

Here, Q{t) = Qicrt) is a function of the random spin configuration eval- 
uated at time t > 0. The left-hand side averages over the perturbed 
dynamics, depending on the hg, and over the initial equilibrium p; the 
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right-hand side averages over the unperturbed dynamics always start- 
ing in p: {Q{t))eq = p{cr) Q{(j) as the equihbrium is time-invariant. 
The hnear correction contains the response function or generahzed sus- 
ceptibihty i?Qy(t, s) which is our object of study. Formally and leaving 
away further decorations, 



h=0 



An interesting case looks at the response in the magnetization itself, 
taking Q = V = ^ a{i) and then 

Rl'iit^s) =/3j2 0<s<t (4) 

is expressible as a space-time-correlation function in the equilibrium 
process. That formula is valid for all times < s < t and for all 
choices of rates that satisfy detailed balance. It is an example of the 
fluctuation-dissipation theorem for finite-time perturbations. The more 
general equilibrium formula of which (jlj) is a special case reads 

Rllit,s) = (3-^jV{s)Q{t))^^, 0<s<t (5) 

which is again true for any choice ([2]) of the rates that satisfies detailed 
balance. A proof of this is easy by applying first-order time-dependent 
perturbation theory and by inserting the equilibrium condition ([1]). 

If we integrate over s G [0, t] with constant hg = h, then 

{Q{t))' - (g(t))eq = hp[{V{t)Q{t)),^ - {V{0)Q{t)U (6) 

Taking t ] +oo we recognize the usual change in the equilibrium 
Boltzmann-factor to first order in h when changing the potential U — >■ 
U-hV. 



The conclusion in equilibrium: adding a potential to a system with 
an equilibrium dynamics is unambiguous, at least when looking at lin- 
ear response and there is a simple and explicit linear response formula 
in terms of an equilibrium correlation in which we recognize the Boltz- 
mann factor. 

The rest of this paper addresses the question what happens if the 
unperturbed dynamics is out-of-equilibrium. The answer is again an 
explicit formula (see (|T2 i) -(l20 l) below) but the choice of how to add 
a potential now does have some influence on the response formula. 
Moreover, there is an interesting interpretation of the resulting corre- 
lations in terms of dynamical fluctuations, which extends equilibrium 
considerations — see Proposition [31 
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3. Going nonequilibrium 

The extension of the previous problem to a nonequihbrium set-up 
has been considered in many papers. We take here the approach of 

[HE]. 

We consider a Markov stochastic dynamics for a finite system. Denote 
the state space by K. We have transition rates W{x,y),x,y G K. We 
do no longer assume that there is a potential, i.e. a function U{x),x e 
K for which y) exp — f/(x) = W {y , x) ex.]) —U {y) . In particular, 
for a stationary distribution p{x),x G K, while 

^[p(x)iy(x,?/) -p(?/)iy(?/,x)] =0, xeK 

y&K 

still, there are nonzero currents of the form p{x) W{x, y)—p{y) Wiij, x) ^ 
for some pairs x ^ y & K. The stationary process (Markov dynamics 
in p) is then no longer time-reversible. We have in mind systems of 
stochastically interacting particles which are driven away from equilib- 
rium; the state x is then the total configuration of particles and the 
transitions are local. An example follows in Section 
Secondly, we also do not need to assume that we start at time t = 
from a stationary distribution. Rather, we have an arbitrary probabil- 
ity distribution p{x),x G K, from which the initial data are drawn and 
then for t > we apply the perturbed dynamics. 

The question is first how to perturb the transition rates W{x, y) — *• 
Wf(x, ?/), by adding an extra potential —htV to the system. For con- 
venience we assume that the perturbation h^^s > 0, is twice differen- 
tiable. Our (physical) assumption here is that the perturbed rates at 
time t > should satisfy 

Wt{x,y) ^ Wix,y) iviy)-vix)] /^x 
Wt{y,x) W{y,x) ^ > 

The inverse temperature j3 signals that the perturbation concerns an 
additional energy exchange with a reservoir at temperature (3~^. The 
assumption ([7]) is conform the condition of local detailed balance as 
often applied in particle systems. That is why we speak of an energy 
impulse. 

Condition (JTj) leaves many possible choices for the perturbed transi- 
tion rates. A quite general choice is 

Wt{x, y) = W{x, y) e'^*[^^(3/)-'^^(-)] (8) 

where the a, 6 G M are independent of the potential V; it was considered 
in [1]. To satisfy ([7]) we need that a + b = P but a or 6 can still vary. 
A first choice is 

Vr/'^ (x, y) = W{x, y) e^[^(s')-^(-)l (9) 
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in which case a = b = f3 /2; that is sometimes called the force-model 
and was explicitly treated in [2]. A second case is 

M/f ) (x, y) = Wix, y) e-^^'^(") (10) 

whence 6 = 0, a = /3, or the opposite 6 = /3, a = 0. 

It is instructive to understand the difference between these cases: we 

can rewrite ([H]) as 

Wt{x,y) = l^(x,y)e'^''^[^(^)+^(^)le^[^(^)-^(^)l (11) 

and we see that making a ^ h gives an extra x ^ y symmetric but 
time-dependent prefactor '0t(a;, y) = ^t{y, = exp ht ^-^{V{y) + V{x)) 
with respect to the force-model of ([9]). Visualizing the situation in 
terms of a one-dimensional potential landscape we imagine the states 
X located at the local minima of a potential U and separated from 
each other via energy barriers. The Arrhenius formula then predicts 
a rate W{x,y) oc exp —P[D{x,y) — U{x)] where D{x,y) = D{y,x) is 
the barrier height between states x and y. Naturally, adding a time- 
dependent potential landscape can affect both the symmetric prefactor 
D{x,y) and the local minima U{x) themselves which gives a possi- 
ble interpretation of the two constants a and b. For example, choos- 
ing ( |TOll only changes the depth of the local minima (binding energies 
U{x) — > U{x) — htV{x)) and not the barrier heights. The above picture 
works best under equilibrium conditions, but one now imagines that 
the nonequilibrium driving adds further asymmetries. 

The linear response question remains unchanged: at time t > the 
expected value of an observable Q will probably deviate from the ex- 
pectation under the unperturbed dynamics. Linear response theory 
out-of-equilibrium is interested in estimating and interpreting the de- 
viations 

to first order in h. We have abbreviated Q(t) = Q{xt) for the observable 
at time t. In other words, we want to compute Rqylt, s) = R{t, s), < 
s < t, in 

= (Qit)), + r ds hs R{t, s) + oih) 

Jo 



4. Response formula 

The present section computes the response function R{t, s) for the 
general perturbation of the form ([H]). In the formula appears the back- 
ward generator L of the jump process; in terms of the transition rates, 



s=0 



Proposition 1. For a perturbation of the form the response func- 
tion is equal to 

R{t, s) = 6^ {Vix,)Q{x,))^ - {V{xs)Qix,))^ 



+ 



{V{x,)LQ{xt)), - {LV{x,)Q{x,))^ (12) 



The proof of this result is essentially a linear order perturbation 
of the Girsanov-formula for the density of the perturbed versus the 
original path-space measures. 

Proof of PropositionUi To see where we must go, we first rewrite the 
right-hand side of ([12]) • In particular, the third term involving LQ can 
directly be combined with the second term, time-derivative in t, adding 
up to (6 — a) multiplied with 

^^{Vixs)Qixt)), = -j-^{Vixs)Qixt)), + Y,f^s{x)Vix)e^'-^^'^Qix) 

X 

where /is(x) = J2y^i.yy^)f^s{y) — J2y^i^^y)f's{.x) solves the master 
equation starting from /xq = A*. As a consequence, we really must prove 
that 

d_ 

ds 

-b{LV{xs)Q{xt))^ - {b-a)Y,{W{xs^y)V{xs)Q{xt))^ (13) 



R{t,s) = a—{V{xs)Q{xt))f, 



+ (b-a) Y,f^s{y)W{y,x)V{x)e'^'-'^''Q{x) 

x,y 



Let a path be denoted by a; = {xs)s, s G [0,t], for Xs G K. Paths are 
piecewise constant and chosen with left limits and right continuous at 
every jump time. For the perturbed process 

r dP^ 

m))l = J dP,{u)^^{u;)Q{x,) (14) 

where we have inserted the density between the path- measures PI^{ijj) 
for the perturbed and the unperturbed {h = 0) Markov dynamics start- 
ing from law fi at time zero. Explicitly (see e.g. Appendix 2 in [6]), 



the Girsanov formula gives 

log-^H = ^ /i46\/(a;,)-ay(a;,-)] 



s<t 

- V fds W{Xs, y) [e^s[bV{y)-aV{xs)] ^^5^ 

y 

where the first sum is over all the jump times s G [0,t]. Up to linear 
order in and with some reordering of the terms, this becomes 

log^M = {h-a)Y,hsV{xs)+aY,hs[V{x,)-V{x,-)] 

^ s<t s<t 

-b [ dshsLV{xs) - (b^a) [ ds /i, V 

Jo Jo y 

(16) 

Higher order in h can easily be controlled. Its second term on the right 
still allows a partial summation into 

J2hs[Vixs)-V{x,-)] = htV{xt)-hoVixo)-Y,yi^s-)[hs-K-] 

s<t s<t 

r' d 

= htVixt)-hoVixo) - / ds—hsVixs) (17) 

Jo ds 

The expression ( IT6l) must now be multiplied with Q{xt) and averaged 
over the original Markov process starting from /i, after which we note 
that 

{{[htV{xt)-hoVixo)]-j ds-^hs Vixs)}Qixt))f. = j dsK^^{VsQt) 

' (18) 

reproduces the first term in fll3p . The last two terms in 0161) are also 
easily identified giving rise to the two middle terms in (fT^ . That leaves 
us with the very first term in fll6p for which must hold that 

{Y,hsV{xs)Q{xt))^. = Y^ I dshsiis{y)W{y,x)V{x)e^'-'^'^Q{x) 

s<t x,y •'^ 

That is indeed true as can be seen by writing the sum over all jump 
times in terms of the random measure dks{y, x) on paths u, which gives 
1 when there is a jump y —>■ x at time s, and is zero otherwise: 

{Y,hsVixs)Q{xt))^ = J2vi^) I K{dK{y,x)Q{xt))^ 

s<t x,y 

By the Markov property e^^~^^^Q{x) = {Q{xt)\xs = x,Xs- = y)^ and 
{dUy, x) Q{xt)), = fiM Wiy, x)e(*-^)^g(x) ds 
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SO that the conclusion fll2p is reached. □ 

5. Example 

We come back to the example ([2]) of a purely dissipative spin-flip 
dynamics. We now add a mixing dynamics. More speciflcally, we not 
only have transitions a with corresponding rates W{a,a^), but 

now we also allow transitions a — > a^^ where the spins at neighboring 
vertices z ~ j' G A get exchanged: a^^{k) = a{k), ii i ^ k ^ j while 
cr*-'(z) = cr'^^^i) = cr{i). The rate for these exchanges is A > 0. The 
result is a reaction- diffusion process on K = {+1, —1}^ with generator 
L acting on functions f : K ^ M., 

Lf{a) = J2 Wia, a^ifia^) - /(a)] + A - f{a)] 

jeA i~j 

That unperturbed dynamics does not satisfy the condition of detailed 
balance when /3 7^ for a nontrivial energy function U (a) in ([2]). There 
is a stationary distribution p of which very little is known; in particular 
it can depend on the -j/^ in ([2]). 

We still consider the magnetization V{(r) = Q{(r) = orga- 
nizing and evaluating the perturbation of amplitude ht,t > 0. Note 
that V{a^^) = V{a) and the transition a — > a^^ leaves the total mag- 
netization unchanged. Hence LV{a) = — 2 (jjiy((T, a*) is still the 
dissipation of magnetization due to the spin flip reaction. Let us ab- 
breviate Ji{cr) = —2a{i) W{a,a^) for the systematic rate of change in 
the local magnetization. We get the linear response around steady 
nonequilibrium from ( |T2l) : 

d 



d 



We see that the equilibrium expression (j4]) gets modifled by the corre- 
lation between at{j) and the flux Ji{as). For a constant perturbation 
hg = h, s ^ [0,t], we can integrate over s G [0,t] to get the leading 
order of the response: 



- 6 V / ds{Ji{ao)asij))p 

JO 



i,jeA- 

Note that the rate A is hiding in the correlation functions but the form 
( ITQI) is unchanged no matter what is A. 

The example is a more microscopic version of a reaction-diffusion model 
but it can also be considered as a toy model for a granular lattice gas 
undergoing inelastic collisions. The spins refer then to the presence 
or absence of energy packets which diffuse but can also be created or 
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get lost. The latter specifies the temperature of the environment and 
the fiux Ji in the above would be the systematic rate of local energy 
change. 

6. More symmetries 

In a number of cases the response formula ( fT2|) simplifies. 
There is first the case where the initial distribution is the stationary 
measure p. Then, i.e., when /i = p is the stationary distribution, 
correlation functions like (y {xs)Q{xt)) p are functions of t — s, so that 
the response function becomes 

RQv{t, s) = a^^{V{xs)Q{xt))p - b{LV{xs)Q{xt))p (20) 

In equilibrium, i.e., under time- reversal symmetry, the two terms in the 
right-hand side of fl2UI) coincide and we recover whenever a + b = /3 
(and independent of ip in (^). 

For the case ([9]), this means that 6 = a = |, the response formula 
becomes 

RQvit,s) = ^^^{Vix,)Qixt))p - ^{LV{xs)Qixt)), 
which is in exact agreement with [1]. 

A special case arises when b = and a = /5 in (jH]), because then the 
response is of the same form as in equilibrium: 

R{t,s) = -/3^^{V{xs)Q{xt)), 

This is indeed a special kind of perturbation, as can also be seen from 
the following consideration. Take h to be constant; the law defined 
by p'^{x) oc p{x)e^'^^^^^ is stationary for the new dynamics (to all orders 
in h) . In other words, here the resulting behavior under this perturba- 
tion is like in equilibrium, even though the unperturbed dynamics can 
be far from equilibrium. 

That last remark brings us to considering the limit t | oo of ( fT2l) in 
which the response formula should show stationary response. Imagine 
thus that we apply a new time-independent dynamics with rates 

W^{x, y) = W{x, y) e''[^^(?')-'^^W], small constant h (21) 

We assume that both the original {h = 0) and the perturbed dy- 
namics show exponential ergodicity in converging to p, respectively 
p^ . Similarly we can replace in the above the function V by another 
function M on and construct p^. Both p^ and p'^ depend on h 
and we investigate their change with respect to the original p to first 
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order in h. The original backward generator is still L. The following 
proposition looks at a special observable, and we write 



Xmv - 



(22) 

The dependence on the constants a, b, h is not made explicit but sits 
in the perturbed dynamics, as in ( I2T1) for perturbing potential V and 
similarly for perturbation M. 

Proposition 2. The stationary response functions ( l22l) equal 

Xtv = XvM = b{MLV), + a{VLM), (23) 

Observe the symmetry when interchanging M and V together with 
the exchange of a and b. If the perturbation is of the form ([9]), then 
only interchanging M and V is enough. This symmetry appears useful 
because it reduces the amount of response functions to be measured. 
Moreover, some experimentally difficult responses can be made more 
accessible by interchanging the role of observable and perturbation. 
While its proof is trivial, we are not aware that this symmetry (1231) has 
been observed before. On the level of generators it simply amounts to 
the direct observation that 

(Lr, - L)M = {Lti - L)V + h{b - a) L{MV) + 0{h') (24) 

where for example L^f^f^x) = W{x, y) eic^[h{bV {y) — aV {x))] [f{y) — 
/(x)]. The symmetry in ( 123|) then easily follows from averaging the 
identity ( 12^ . We add however a different proof that connects with the 
response formula ( fT2l) . 

Proof of Proposition [B For observables Q which are of the form Q = 
LM, the linear response is given by (fT2|) . 



^LM,y(^'^) = §-Jci{V{xs)LM{xt))p + b{LVixs)M{xt))p] (25) 

When we consider a constant perturbation hg = h, we get the inte- 
grated form of the response function: 

R{t) = f <lsR%,y{t,s) 
Jo 

= a{[Vixt) - Vixo)]LMixt)), + b{[LV{xt) - LV{xo)]Mixt)), 



It suffices to take t ^ cxd to see the appearance of ( 1231) . The exchange 
of the limits h ^ and t f +oo is trivial in the case considered, so 
that for constant ht = h 



hmR{t) = xtv 



□ 



12 



7. Response and dynamical fluctuations 

We turn to the interpretation of the response functions in terms of 
fluctuation theory. The standard interpretation of the equihbrium re- 
sponse ([5]) is in terms of energy dissipation, and that is why ([5]) is called 
the fluctuation-dissipation theorem even though it really deals with re- 
sponse. That terminology and corresponding interpretation remains 
true and useful for the first term in ( 120|) at least when considering the 
flux in excess to what already was present (since we now deal with 
nonequilibrium) . That was explained in [2], section 5, and we also see 
it in the identity ffTBj) which can be interpreted as a conservation of 
energy. From a probabilistic point of view it is more interesting to con- 
centrate on what is new with respect to equilibrium, the second term 



The way of responding and the way of fluctuating are like each other's 
time-reversals. For inspiration, we turn again to the equilibrium (hence, 
time-reversible) case, where the response to a perturbation typically 
goes along the same path as that of a spontaneous fluctuation; that is 
sometimes called Onsager's regression hypothesis and in our context it 
could be summarized as L = L* where L* is the adjoint in the p-scalar 
product. In nonequilibrium the regression of a fluctuation is also the 
time-reversal of its appearance, but now the time-reversal is not trivial. 
In particular the second term in ( 120|) is 



which now cannot be written as a time-derivative as in the first term 
of (pUj) . Yet, it is related to a fluctuation, as we explain now. 
Suppose a (constant) perturbation V is added to the system. The 
system responds but in the long time, the perturbation also installs a 
new stationary law. Let us denote this new stationary law by /i. On the 
other hand one can compute the probability that in the unperturbed 
dynamics ^ occurs as a fluctuation. This takes us to the dynamical 
fluctuation theory for Markov processes, started by [5], see also e.g. 
in [3]. Without going to the full details it suffices here to recall that 
for an ergodic Markov process with backward generator L there is a 
fluctuation functional /(/x) on the probability laws /i on of the form 



in ([20]). 



h ^p{x)V{x) L*e(*-^)^g(x) 



X 




(26) 



for which in the sense of the theory of large deviations 



Vmh p[pr ~ /i] ^ e T t +00 
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for the empirical distribution 

1 r 

Pt{x) = — 6xt,x dt, with 6a,b = if a 7^ fe and Sa,b = 1 if a = fe 



of occupation times over the time-interval [0,r]. We refer to [9] for 
a dynamical fluctuation theory in the context of the present paper. 
In a sense, exp — r/(/i) gives the plausibility of the long-term (= r) 
appearance of (= dynamical fluctuation to) the statistics /i. Taking 
g = exp{bhM/2) in ([211), we see 



J2 W{x, y)e'tmy)-Mix)] _ J2 W{x, y) 



(27) 

That M can now be interpreted as a potential. The infimum in flTTj) 
gets reached a.t M = V, the potential for which fi is the stationary 
law. Already here we see a complementarity between response and 
fluctuations: a perturbation V gives a new stationary law /i, and to 
flnd the probability of a fluctuation fi in the original dynamics, one 
has to flnd exactly this V. But there is also a quantitative relation, 
in particular as realized in the second term in (120|) . as we prove in the 
next proposition. 

Proposition 3. The dynamical fluctuation functional I[ix) satisfies 

X 

Proof of Proposition 0. For our case we can exchange the limit by the 
inflmum and the small h limit, see also [8]. We can then compute (l26l) 
to flrst order in h by taking M = bhV/2 in (1271) . and expanding 

J2fi{x) W{x,y) [1 - e^Mny)-v^W]/2] ^ ^;.(x) W{x,y)[V{y) - V{x)] 

x,y x,y 

x,y 

where we have used already that {LV)p = by stationarity. For the 
second term we can replace the p(x) by fi{x) because we are already at 
second order in h and write 

E /^(^) ^(^' ^) (^(^) - ^(^))' + ^(^') = 

x,y 



^ 5^/.(x) W{x,y)[Viy) - V(x)] [1 - e^Mn.)-^^)]/^] 



6/1 

^-^J2^^(x)Wix,y)[Viy)-Vix)] 

x,y 
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because n is invariant under the perturbed dynamics. Collecting all 
terms we get the result (12^ . 

□ 

Of course the LV{x) in the response formulae ( fT2l) -(l20 l) -( l23|) has the 
usual meaning of being the expected rate of change in V while at x. 
Proposition [3] adds the interpretation that it can also be seen as the 
change in escape rate from x when adding a potential V. From fl28p 
the second term in fl20l) gives a correlation with a generalized escape 
rate and thus relates with the dynamical fluctuations of the occupation 
times, PE]. 



8. Conclusion 

We have generalized the results of P, [2] to the perturbation first con- 
sidered by [1], Proposition [H We have also added a stationary response 
relation and noted a new symmetry, Proposition [21 The fluctuation in- 
terpretation of [2] remains intact. Proposition [31 
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